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The notion of a Fermi surface (FS) is one of the most ingenious concepts 
developed by sohd state physicists during the past century. It plays a central 
role in our understanding of interacting electron systems. Extraordinary efforts 
have been undertaken, both by experiment and by theory, to reveal the FS of 
the high temperature superconductors (HTSC), the most prominent strongly 
correlated superconductors. Here, we discuss some of the prevalent methods 
used to determine the FS and show that they lead generally to erroneous results 
close to half filling and at low temperatures, due to the large superconducting 
gap (pseudogap) below (above) the superconducting transition temperature. 
Our findings provide a perspective on the interplay between strong correlations 
and superconductivity and highlight the importance of strong coupling theories 
for the characterization as well as the determination of the underlying FS in 
ARPES experiments. 

During the last decade, Angle Resolved Photoemission Spectroscopy (ARPES) has 
emerged powerful tool to study the electronic structure of the HTSC This is 
because ARPES is a direct method to probe the FS, the locus in momentum space where 
the one electron excitations are gapless [4] . However, since the low temperature phase of the 
HTSC has a superconducting or pseudogap with d-wave symmetry, an FS can be defined 
only along the nodal directions or along the so-called Fermi arcs 
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. The full 'un- 
derlying FS' emerges only when the pairing interactions are turned off, either by a Gedanken 
experiment, or by raising the temperature. Its experimental determination presents a great 
challenge since ARPES is more accurate at lower temperatures. Since the FS plays a key 
role in our understanding of condensed matter, it is of importance to know what is exactly 
measured by ARPES in a superconducting or in a pseudogap state. The problem becomes 



even more acute in HTSC due to the presence of strong correlation effects j8|, |9|, |10|, 



Hence, it is desirable to examine a reference d-wave superconducting state with aspects of 



w these considerations, we 
111 and discuss our results 



strong correlation built explicitly in its construction. Motivated 
study the FS of a strongly correlated d-wave superconductor p 
in the context of ARPES in HTSC. 

We begin by highlighting the differences between a Fermi and a Luttinger surface. The 
FS is determined by the poles of the one electron Green's function The Luttinger 
surface is defined as the locus of points in reciprocal space, where the one particle Green's 
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function changes sign [12] . In the Fermi hquid state of normal metals, the Luttinger surface 
coincides with the FS. In a Mott-Hubbard insulator the Green's function changes sign due 
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at momenta k of 



to a characteristic 1 /cu-divergence of the single particle self energy 
the non-interacting Fermi surface. In the HTSC the gapped states destroy the FS but only 
mask the Luttinger surface. Hence, it seems natural to relate the Luttinger surface of the 
superconducting and of the pseudogap states with the concept of an 'underlying FS', and 
ask if such a surface can be determined by ARPES. 

To answer this question, we recall that the elementary excitations in a superconductor 
are given by the dispersion relation. 



Ek = -^el + Al, ^, = ek-fi, (1) 

where are the momentum dependent orbital energies of electrons in the absence of a super- 
conducting order parameter Aj^; fi is the chemical potential. The corresponding Luttinger 
surface is determined by the condition C,k = 0, which is also the definition of the normal state 
FS when = 0. In the following, we discuss two methods commonly used to determine 
the underlying FS, viz., the Luttinger surface, of the HTSC by ARPES jll B ll^. 

In the so-called 'maximal intensity method' the intensity of ARPES spectra at zero 
frequency is used to map out the underlying FS. It can be shown that this quantity is, 

^ k 

where, F^ is determined both by the experimental resolution and the width of the quasipar- 
ticle peak. When the momentum dependence of F^ is small compared to that of (as is 
usually the case), the maximal intensity is given by the set of momenta hk for which E^ is 
minimal. 

To examine the accuracy of this method in determining the underlying FS, we calculate 
this quantity for a strongly correlated (i-wave superconducting state. All calculations are 
done with model parameters for HTSC using the renormalized mean field theory (RMFT) 



111 ], for which the quasiparticle dispersion E^ retains the form of Eq.^ In Fig- 1, we show 
our results for the spectral intensity at zero frequency as well as the locus of the Luttinger 
surface. The former is deduced from the inverse of Ek- 

For large hole doping, x = 0.25, the superconducting gap is small and the Luttinger 
surface is close to the points in momentum space for which the zero frequency intensity is 



maximal. But for smaller doping, x = 0.05, the gap is substantial and the Luttinger surface 
deviates qualitatively from the maximal intensity surface due to the momentum dependence 
of Afc (see ridges in Fig. 1). We have verified that this behaviour persists for a wide range 
of |Afc|, and not just the values estimated from RMFT used in Fig. 1. Although not widely 
discussed in the literature, this splitting may be deduced from experimental data, e.g., the 
intensity plots in E — k space along symmetric lines (0, 0) — (vr, 0) — (vr, vr) in It follows 
that when the gap or the pseudogap is large, the criterion of maximal spectral intensity alone 
does not suffice to identify the correct FS and it is necessary to supplement the analysis of 
the zero frequency ARPES intensity (Eq. 2) with a dispersion relation such as Eq. 1. These 
considerations explain why the (outer) maximal intensity ridges seen in ARPES (at low 
temperatures in the underdoped regime) may yield an underlying FS whose volume is too 



large 



Another method used in extracting the Luttinger surface is the 'maximal gradient 
method'. The method is based on the fact that the FS is given by the set of /c- values 
for which the momentum distribution function Uk shows a jump discontinuity. When this 
discontinuity is smeared out, say, by thermal broadening or a small gap, the gradient of Uk, 
|Vnfc|, is assumed to be maximal at the locus of the underlying FS. 

We calculated \Vnk\ within RMFT and show our results in Fig. 2. We see that the 
maximal gradient surface is very sensitive to the presence of even small gaps. For example, 
the superconducting gap at x = 0.25 is quite small. Nonetheless, the electron-like Luttinger 
surface (determined by ^k = 0) is not clearly revealed by the ridges in |Vnfc|. Similar 
deviations of |Vnfc| from the underlying surface are also obtained from a high temperature 
expansion of the t — J model and dynamical cluster approximation in the Hubbard model 
19| . We conclude that the maximal gradient method alone cannot be used to determine the 
underlying FS unambiguously from ARPES data. 

The notion that the underlying FS of a pseudogapped or a superconducting state is 



identical to the Luttinger surface is only approximately correct |12l. l20l|. In the Fermi liquid 
state of normal metals, the FS satisfies the Luttinger sum rule; the volume enclosed by the 
FS is identical to the total number of conducting electrons. But, in a superconductor, the 
chemical potential is generally renormalized and is a function of the superconducting order 
parameter, fi = /^^^(A). The number of states niutti^) enclosed by the resulting Luttinger 
surface, C,k = 0, then deviates from the true particle number n, as the results in Fig. 3 show. 
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However, this effect is small (a few percent) and unlikely to be discerned experimentally. 
The discrepancy between riLutti^) and n vanishes when particle-hole sjTiimetry is present. 
Further, it changes sign when the geometry of the Luttinger surface changes from hole-like 
to electron-like, as seen in Fig. 3. 

Finally, we focus on the influence of the strong electron-electron interactions on the 
geometry of the Luttinger surface close to half filling. The Cu-0 planes of the HTSC are 
characterized by a nearest neighbor (NN) hopping parameter t ~ 300 meV and a next nearest 
neighbor (NNN) hopping parameter t' ^ —t/4. These parameters are the bare parameters, 
and determine the dispersion relation, 

efc = —2t(^cos{k.j.) + cos{ky)^ — 2t' (^cos{kx + ky) + cos{kx ~ ky)j , (3) 

in the absence of any electron-electron interaction. On the other hand, true hopping pro- 
cesses are influenced by the Coulomb interaction U ^ 12t leading to a renormalization of 
the effective hopping matrix elements, 

t i=i{U), t' ^ i = i'{U). 

Close to half filling we find t cc J = At'^/U and t' —>■ 0, i.e., the frustrating NNN hopping is 
renormalized to zero. This behavior is illustrated in Fig. 4. The resulting Luttinger surface 
renormalizes to perfect nesting. A similar behavior has been observed in recent variational 
studies of organic charge transfer-salt superconductors [21]. At half filling, the so-called 
Marshall sign rule is valid rigorously in the absence of frustration. The degree of effective 
frustration can then be estimated by the size of deviation from the Marshall sign rule as a 
function of the (frustrating) bare t'. A numerical study has found, that the Marshall sign 
rule remains valid even for small but finite t', viz., the effective frustration renormalizes to 
zero Q]. This behavior is in agreement with the results presented in Fig. 4, and is unique 
to strong coupling theories such as RMFT. 

We showed that the accurate determination of the underlying FS in underdoped HTSC 
is a difficult task and that analysis of the experimental data alone is often insufficient for 
an unambiguous determination of the FS. Commonly used methods like the zero frequency 
spectral intensity or the gradient of can yield significant deviations from the true Luttinger 
surface as shown in Fig. 1 and 2. Indeed, a clear distinction between electron- and hole-like 
underlying FS cannot be made solely from analyses of spectral intensity maps when the gaps 
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are large. Such analyses have to be supplemented by a minimal modelling of the gapped 
states. Furthermore, the underlying FS in the pseudogapped or superconducting state fulfills 
Luttinger theorem only approximately, owing to the dependence of the chemical potential on 
the superconducting gap. We also demonstrated that the strong correlations renormalize the 
ratio t' / 1 near half filling, yielding a Luttinger surface which is perfectly nested. This suggests 
in a very natural way that the strong coupling mean field superconducting state is unstable 
to antiferromagnetism at low doping. Our findings resulting from the combined effects of 
strong correlations and d-wave superconductivity, allow for a more precise interpretation of 
experiments that determine the FS of HTSC 
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FIG. 1: The zero frequency spectral intensity (deduced from the inverse of E^) in the first Brillouin 
zone for hole dopings x = 0.05 (left) and x = 0.25 (right). The color coding blue/red corresponds 
to the low/high zero frequency spectral intensity. The ridges of maximal intensity are indicated by 
the (dashed) red and (dashed-doted) orange lines respectively, the Luttinger surface by the black 
line. 




FIG. 2: The gradient of the momentum distribution function, \nk\, in the first Brillouin zone, for 
hole dopings x = 0.05 (left) and x = 0.25 (right). The color coding is blue/red for small/large values 
of |Vfcnfc|. The ridges of maximal |Vjtnfc| are indicated by the (dashed) red and (dashed-doted) 
orange lines respectively, the Luttinger surface by the black line. 
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FIG. 3: The deviation, riLutt — n, of the actual volume of the Luttinger surface from the Luttinger 
sum-rule, as a function of hole-doping x. Calculations are performed by RMFT {t' = —t/A, 
U = 12t). The deviation is minimal when the topology of the Luttinger surface changes from 
hole-like to electron-like. 

Inserts (a) and (b): Model calculation for the renormalization of the chemical potential, fJ-sc ~ fJ-FS: 
and the resulting nLutt — n as a function of the d-wave order A, for various doping x. 



9 




(0,71) 



(7i;,7i;) 



0.05 
hole doping x 



V >' 1 ' 1 ' 1 
" Vv (b) 

— \v\ 

\\\ 
\\\ 
\\ \ 
\ \^ \ 

— v\\ \ 


— x=0 
-- x=0.08 

— x=0.16 
-- x=0.24 

— x=0.32 


\\\\ \ 

\ 

\>\s \ 

— nSn.^ \ 
\s\n 

\s\ 

(0,0) 

\ 1 , 1 , 1 


\ 

, 1 



FIG. 4: (a) Renormalization of the next nearest neighbor hopping ampHtude, t' — > t', as a function 
of hole doping concentration x for various values of bare t' . All effective t' are renormalized to zero 
at half filling by the large Coulomb repulsion. We highlight the region for which we expect the 
superconducting d-wave state to become unstable against antiferromagnetism (AFM) due to the 
nearly perfect nesting of the Luttinger surface. 

(b) The geometry of the Luttinger surface for the high temperature superconductors {t' = —t/4). 
The change is non-monotonic for small doping x, when the Luttinger surface is renormalized to 
perfect nesting due to the strong Coulomb interaction. For x = 0.16, the topology of the Luttinger 
surface changes from hole-like to electron-like. 

Calculations are performed for the Hubbard model with U = 12t, using RMFT. 
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